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ABSTRACT 

Context. Continued investigation of the linkage between magneto-acoustic energy generation in stellar convective zones 
and the energy dissipation and radiative emission in outer stellar atmospheres in stars of different activity levels. 
Aims. We compute the wave energy fluxes carried by longitudinal tube waves along vertically oriented thin magnetic 
fluxes tubes embedded in the atmospheres of theoretical main-sequence stars based on stellar parameters deduced by 
R. L. Kurucz and D. F. Gray. Additionally, we present a fitting formula for the wave energy flux based on the governing 
stellar and magnetic parameters. 

Methods. A modified theory of turbulence generation based on the mixing-length concept is combined with the magneto- 
hydrodynamic equations to numerically account for the wave energies generated at the base of magnetic flux tubes. 
Results. The results indicate a stiff dependence of the generated wave energy on the stellar and magnetic parameters 
in principal agreement with previous studies. The wave energy flux Jltw decreases by about a factor of 1.7 between 
GO V and K0 V stars, but drops by almost two orders of magnitude between K0 V and M0 V stars. In addition, the 
values for Fltw are significantly higher for lower in-tube magnetic field strengths. Both results are consistent with the 
findings from previous studies. 

Conclusions. Our study will add to the description of magnetic energy generation in late-type main-sequence stars. Our 
results will be helpful for calculating theoretical atmospheric models for stars of different levels of magnetic activity. 

Key words, methods: numerical — MHD — stars: chromosphere — stars: photosphere — stars: magnetic fields — 
waves 

1. Introduction K2 V stars with different levels of magnetic activity with 

■ the filling factor for the magnetic component determined 

G> An outstanding problem in stellar astrophysics concerns from an observational relationship between the measured 

O i the identification of physical processes responsible for the magnetic area cover age and the stellar rotation period. For 

O . heating of outer stellar atmospheres and the acceleration starg with yery glow rotatiorlj they were able to repr oduce 

^ehar wind s (see review s by | Naram fc Ulmschneideij the bagal flux Hmit of chromospheric emiss ion previously 

[19901 119961 and | Gudel || 2007| ). For the Sun and other types identi fied with no n -magn etic regions. Most notably, how- 

of stars with surface convection zones, acoustic heating ^ | Cuntz et aL | (f^ de duced a relationship between 

has been identified as most likely responsible for balm- the Ca n H+K emission and the stellar rotat i n rate that 

mg the basal flux emission (e.g., | Buchholz et al. | [1998|; ig consistent with th e relationsh ip previously obtained by 

i Cuntz et al. || 2007| ). On the other hand, it is well known that observa ti ns; see also lCuntz et a l. (1998) for earlier results, 
most, if not all stars also exhibit a large amount of magnetic 

activity. Thus, the chromospheres of main-sequence stars, Further studies for a large spectral range of stars were 

including the Sun, are expected to be significant ly shaped S lven by |Fawz Y et al. | ((2002) based on specified values for 

bymagnetically heated structure (e.g. Jaar 1994; Sch^S the magnetic filling factor. They concluded that heating by 

199G ) acoustic and longitudinal flux tube waves is able to explain 

There is a large body of previous work devoted to the most of , the , observed range of chromospheric activity as 

description of the two-component structure of stellar chro- S a ^ ed h ? \ he Ca 11 and M § U , lmes - 0n the other hand ' 

mospheres. In these models the magnetic component of mdlrect evidence for non-wave (i.e., reconnective) heating 

the chromosphere is typically he ated by energ y dissipation ™* also deduced needed to exp am the structure of the 

of longitudinal flux tube waves. ICuntz et alT (iSh com- hl § hest layers of stellar chromospheres, 
puted two-component theoretical chromosphere models for This type of models, as well as envisioned future models 

of chromospheric heating and emission, partially motivated 

Send offprint requests to: Diaa E. Fawzy by the quest of investigating the effects of UV and EUV 
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emission on planetary atmosphere and (potentially) the 
evolution of life (e.g.. iGuinan et al. I 120031 : iLammer et al.l 
120031: iGiidel 1 120071: ICuntz et alJl2010D . require the continu- 
ation of detailed simulations of magnetic wave energy gen- 
eration, including studies on longitudinal tube waves in dif- 
ferent types of stars, particularly main-sequence stars. This 
latter goal is the focus of the present paper. 

Previous work on the calculation of lo ngitudinal tube 
waves has been based on progress made by iMusielak et ail 
(|l994f) who corrected the Lighthill-Stein theory by incor- 
porating an improved description of the spatial and tem- 
poral spectrum of the turbulent convection and utilized 
the corrected theory for calculating revised stellar acous- 
tic wave energy fluxes (|Ulmschneider et all Il996l [l999) . 
This type of work focused on the generation of acous- 
tic waves; however, it did not consider stellar magnetic 
fields. Considering the fundamental importance of mag- 
netic heating in most, if not all stars, a set of papers fo- 
cused on the study of longitudinal and transverse tube 
wave generation has been pursued (e.g. IMusielak et al.l 
119891 119951: lUlmschneider fc Musielakll 19981 ). In subsequent 
work, lUlmschneider et al.l (1200 if) u s ed the approach devel- 
oped by lUlmschneider fc MusielakI (|1998t ) to compute the 
wave energy fluxes carried by longitudinal tube waves prop- 
agating along thin and vertically oriented magnetic flux 
tubes that are embedded in atmospheres of late-type stars. 
This numerical appro ach supplemented previous work by 
IMusielak et al.l (|2000| ). who analytically calculated the lon- 
gitudinal wave energy fluxes generated in stellar convective 
zones. 

In the numerical approach bv lUlmschneider fc MusielakI 
(1998), longitudinal tube waves are generated as a result of 
the squeezing of a thin, vertically oriented magnetic flux 
tube by external pressure fluctuations produced by the 
turbulent motions in a stellar photosphere and convection 
zone, which correspond to associated velocity fluctuations. 
Hence, to compute the pressure fluctuations imposed on the 
tube, it is required to know the external turbulent motions. 
The motions are modeled by specifying the rms velocity 
amplitude and using an extended Kolmogorov turbulent 
ener gy spectrum with a modified Gaussian frequency fac- 
tor (jMusielak et al.lll994f ). 

The main advantage of this approach is that it is not 
restricted to linear waves and that it allows for occasion- 
ally large-amplitude wa ves observed on the Sun at the 
photospher i c lev e l (e.g., Muller _ 19851 : iKomm et al.l 119911 : 
ikesis et al.l 119931: iMuller et al.lll994l) and also seen in de- 
tailed time-d ependent simulations of so l ar and stellar con- 
vecti o n (e.g., Nordlund k, Dravinsl 19901: Nordlund fc Stein 



1991 



1997) 



ICattaneo et al.l Il99lb ISteffenl 119931: iNordlund et al 
Horizontal flow patterns are a notable candidate 
process for the initiation of wave modes wit h respect to 
flux tubes (see Fig. 1); see lStein et al.l (|2009al ibh for recent 
models of convective flows for the Sun based on up-to-date 
simulations extendin g toward the scale of supergranul es. 

The code used bv lUlmschneider fc MusielakI ( 1998h was 
originally developed by iHerbold et al.l (|1985f) who treated 
magnetic flux tubes in the so-called thin flux tube approx- 
imation and described them mathematically by using a set 
of one-dimensional, time-dependent and nonlinear MHD 
equations. It allows to compute the instantaneous and time- 
averaged longitudinal tube wave energy fluxes and the cor- 
responding wave energy spectra. It requires specifying the 
strength of the magnetic field inside the flux tube and the 
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Fig. 1. Diagram of a flux tube embedded into a stellar con- 
vective zone. The squeezing point of the tube is assumed to 
be located at optical depth T5000 = 1> coinciding with the 
"stellar surface". Credit: P. Ulmschneider. 



height in the stellar atmosphere where the squeezing is as- 
sumed to take place. The code has previously been used 
to calculate wave energy fluxes and spectra for longitudi- 
nal tube waves p ropagating in the solar atmosphere (see 
iFawzv et al.lll998l for models of different spreading factors) , 
and to investigate the dependence of these fluxes on the 
magnetic field strength, the rms velocity amplitude of tur- 
bulent motions, and the location of the squeezing in the at- 
mosphere; for r ecent models f or other stars with non-solar 
metallicities see IFawzv! (I2010D . 

The reason for reinvestigating the generation of longitu- 
dinal flux tube waves in main-sequence stars is three-fold. 
First, we would like to use realistic combinations of (T e ff, 
logg), with T c ff as stellar effective temperature and logg 
as surface gravity, for main-sequence stars guided by up- 
to-date studies by R. L. Kurucz and D. F. Gray. Note that 
log o is typically close t o 4.5 (s ee Table 1). Previous models 
by lUlmschneider et al.l (|200ll) and others have been pur- 
sued for either log g — 4 or 5, thus resulting in unnecessary 
interpolation errors. Secondly, we would like to investigate 
the amount of upward propagating longitudinal wave en- 
ergy flux for a wider range of stellar convective and mag- 
netic parameters, notably the mixing length a amid recent 
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Table 1. Theoretical main-sequence stars 



Sp. Type 


T c ff 


R 


M 


log 9 




(K) 


(R@) 


(M ) 




FO V 


7178 


1.620 


1.600 


4.223 


Fl V 


7042 


1.541 


1.560 


4.255 


F2 V 


6909 


1.480 


1.520 


4.279 


F3 V 


6780 


1.453 


1.480 


4.283 


F4 V 


6653 


1.427 


1.440 


4.287 


F5 V 


6528 


1.400 


1.400 


4.292 


F6 V 


6403 


1.333 


1.330 


4.312 


F7 V 


6280 


1.267 


1.260 


4.333 


F8 V 


6160 


1.200 


1.190 


4.355 


F9 V 


6047 


1.155 


1.120 


4.362 


GO V 


5943 


1.120 


1.050 


4.360 


Gl V 


5872 


1.100 


1.022 


4.364 


G2 V 


5811 


1.080 


0.994 


4.368 


G3 V 


5760 


1.037 


0.967 


4.392 


G4 V 


5708 


0.993 


0.940 


4.417 


G5 V 


5657 


0.950 


0.914 


4.443 


G6 V 


5603 


0.937 


0.888 


4.443 


G7 V 


5546 


0.923 


0.863 


4.443 


G8 V 


5486 


0.910 


0.838 


4.443 


G9 V 


5388 


0.870 


0.814 


4.469 


KO V 


5282 


0.830 


0.790 


4.497 


Kl V 


5169 


0.790 


0.766 


4.527 


K2 V 


5055 


0.750 


0.742 


4.558 


K3 V 


4973 


0.730 


0.718 


4.567 


K4 V 


4730 


0.685 


0.694 


4.608 


K5 V 


4487 


0.640 


0.670 


4.651 


K6 V 


4294 


0.601 


0.643 


4.689 


K7 V 


4133 


0.565 


0.614 


4.722 


K8 V 


4006 


0.533 


0.582 


4.749 


K9 V 


3911 


0.505 


0.547 


4.770 


MO V 


3850 


0.480 


0.510 


4.783 



progress made through models bv lStein et aU (|2009al lhh and 
others. Thirdly, we would like to deduce a fitting formula 
for the wave energy flux that allows insight into the role of 
the relevant parameters concerning that flux and, further- 
more, offers a more universal use. 

Our paper is structured as follows: In Sect. 2, we com- 
ment on the parameters of theoretical main-sequence stars. 
Additionally, we summarize the method for the computa- 
tion of longitudinal tube waves as well as the construction 
of stellar flux tube models. Our results are given in Sect. 3. 
Finally, in Sect. 4 we present the summary and conclusions. 

2. Methods 

2.1. Comments on the theoretical main-sequence stars 

Stellar parameter s for t h eoret ical main-sequence stars have 
been deduced by iGravi (|2005l ); see his Table B.l. His val- 
ues, notably T e g and log g, serve as basis for the present 
study. We also improved the accuracy of the log g values if 
more accurate values for the stellar masses and stellar radii 
were given. For stellar spectral types with no data given, we 
calculated those data using biparabolic interpolation. The 
stellar data are summarized in Tabic 1. 

Another set of spectral models has been constructed by 
R. L. Kurucz and collaborators. They take into account 
millions or hundred of millions of lines for a large array 
of atoms and molecules; see, e.g., ICastelli fc Kuruczl ()2004h 
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Fig. 2. Root mean square turbulent velocities at the 
squeezing point (T5000 = 1) for the set of theoretical main- 
sequence stars for different values of the mixing-length pa- 
rameter a. 
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Fig. 3. Power spectra of the instantaneous wave energy 
flux in flux tubes with B/B cq = 0.85 for F5 V, G5 V, 
K5 V, and M0 V stars (from top to bottom) as a function 
of circular frequency uj. The mixing- length parameter is 
assumed as a = 2. 

and|K urucd (|2005h for technical details. These models in- 
dicate ver y simil a r effec tive temperatures compared to the 
models by IGravi (2005) for most types of stars. However, 
stellar spectral types of K5 V and below, the indicated ef- 
fective temperatures of R. L. Kurucz are consistently lower 
noting that the difference amounts to nearly 300 K for spec- 
tral type M0 V. Therefore, we assumed average values be- 
tween the models by D. F. Gray and R. L. Kurucz for stars 
of spectral spectral type K5 V and M0 V in the following. 

2.2. Convective zone models and turbulent velocities 

The method for calculating wave energy fluxes carried 
by longitudinal tube waves ad opted in the presen t pape r 
has been described in detail bv lUlmschneider et~a l. (2001). 
Thus, it is not necessary to present an intricate discussion 
in the following. In the solar application, it is possible to 
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Table 2. Wave energy flux a for different parameters a and 
V 



Sn Tvne 


T 


fur 


-Fltw 


-Fltw 




(k) 












■n = 0.75 


■n = 0.85 


■n = 0.95 


OL 


1.5 








1 V 


6528 


6.03E8 


4.05E8 


9.41E7 


F8 V 


6160 


O . O cfi-JO 


3.10E8 


9.14E7 


GO V 


5943 


4.75E8 


2.94E8 


7.91E7 


G2 V 


5811 


4.50E8 


2.62E8 


7.91E7 


G5 V 


5657 


3.56E8 


2.20E8 


6.53E7 


G8 V 


5486 


3.27E8 


1.98E8 


5.85E7 


KO V 


5282 


2.54E8 


1.56E8 


4.77E7 


K2 V 


5055 


1.92E8 


1.10E8 


3.42E7 


K5 V 


4487 


1.14E8 


6.83E7 


2.12E7 


K8 V 


4006 


4.47E6 


2.87E6 


9.82E5 


MO V 


3850 


2.26E6 


1.33E6 


5.36E5 


OL 


1.8 








1 V 


6528 


8.99E8 


5.00E8 


1.25E8 


F8 V 


6160 


7.51E8 


4.62E8 


1.15E8 


GO V 


5943 


5.69E8 


3.40E8 


9.73E7 


G2 V 


5811 


5.63E8 


3.02E8 


9.05E7 


G5 V 


5657 


5.29E8 


2.78E8 


7.96E7 


G8 V 


5486 


3.98E8 


2.37E8 


7.32E7 


KO V 


5282 


2.93E8 


1.82E8 


5.60E7 


K2 V 


5055 


2.69E8 


1.63E8 


4.72E7 


K5 V 


4487 


1.30E8 


8.15E7 


2.50E7 


K8 V 


4006 


5.36E6 


3.50E6 


1.24E6 


MO V 


3850 


2.28E6 


1.44E6 


6.49E5 


a — 


2.0 








F5 V 


6528 


1.12E9 


6.07E8 


1.47E8 


F8 V 


6160 


8.51E8 


4.92E8 


1.34E8 


GO V 


5943 


6.77E8 


3.96E8 


1.11E8 


G2 V 


5811 


6.35E8 


3.65E8 


1.06E8 


G5 V 


5657 


5.81E8 


3.35E8 


1.01E8 


G8 V 


5486 


4.79E8 


2.86E8 


8.05E7 


KO V 


5282 


3.39E8 


1.95E8 


6.39E7 


K2 V 


5055 


2.98E8 


1.86E8 


5.81E7 


K5 V 


4487 


1.44E8 


9.19E7 


2.86E7 


K8 V 


4006 


6.79E6 


4.37E6 


1.57E6 


MO V 


3850 


3.32E6 


2.11E6 


7.62E5 



"The unit of -Fltw is erg cm s 



select many model parameters and characteristic values di- 
rectly from observations. However, for stars other than the 
Sun such data are mostly unavailable. Therefore, we need 
to discuss in some detail the physical reasoning behind our 
choice of relevant parameters used in our calculations. 

In the current approach, the magnetic flux tubes are 
embedded in nonmagnetized photospheric convection zones 
(see Fig. 1). Considering that the interaction between the 
flux tubes and the convective turbulence is the driving 
mechanism for the generation of longitudinal tube waves, 
among other waves, models of the stellar convection zones 
are required. Guided by previous studies, it is assumed that 
the squeezing of the tube is symmetric with respect to the 
tube axis. The computed pressure fluctuations are subse- 
quently translated into gas pressure and magnetic field fluc- 
tuations inside the tube assuming horizontal pressure bal- 



ance. Finally, the internal velocity perturbation resulting 
from the internal pressure fluctuation is calculated. This 
internal velocity served as a boundary condition in the 
numerical simulation of the generation of the longitudinal 
tube waves. 

Both numerical simulations of stellar convection and 
mixing length models show that the maximum convective 
velocities occur at optic al depths of T5000 ~ 10 to 100. For 
example, ISteffenl (|l993l) found in his time-dependent solar 
numerical convection calculations that maximum convec- 
tive velocities uciiax — 2.8 km s _1 are reached at T5000 ~ 50 
and that these values can be reproduced via mixing length 
theory with a mixing length parameter of a ~ 2. The value 
a = 2 is furthermore indicated by time-dependent hydrody- 
namic si mulations of stellar conve ction for stars other than 
the Sun ()Trampedach et aT1ll997l) as well as by a careful fit- 
ting of evolutionary tracks of the Sun with its present lumi- 
nosit y, effective temperature and age (jSchroder fc Eggletonl 
Il996h . 

Nevertheless, there i s still some debate about th e most 
appropriate value of a. Nordl und fc Dravins! ()1990[) origi- 
nally pursued detailed numerical simulations based on 3-D 
hydrodynamics coupled with 3-D non-grey radiative trans- 
fer for stars similar to Procyon (F5 IV- V), a Cen A (G2 V), 
(3 Hyi (G2 IV), and a Cen B (Kl V). They concluded a 
mixing- length parameter of a = 1.5 (or slightly higher), 
although the mixing length concept appeared to be prob- 
lematic at ph otospheric he i ghts. A mixing length of 1.5 was 
also used by iCuntz et al.l (|1999t ) in their two-component 
theoretical chromosphere models for K2 V stars with dif- 
ferent levels of magnetic activity. Even though the deduced 
relationship between the Ca II H+K emission and the stel- 
lar rotation rate was found to be largely consistent with the 
observed relationship, the agreement could probably be im- 
proved if a somewhat higher longitudinal wave energy flux, 
corresponding to a sligh tly larger m i xing len gth parameter, 
was adopted. Recently, IStein et aD (2009a. h) pursued up- 
dated state-of-the-art simulations of solar convection zone 
extending toward the scale of supergranules indicating a 
mixing length parameter of a ~ 1.8. For these reasons, we 
will calculate a set of models concerning wave energy gener- 
ation of longitudinal tube wave for a set of a values, which 
are a = 1.5, 1.8, and 2.0. 

2.3. Computation of stellar magnetic flux tube models 

Our treatment of stellar convection associated with the fa- 
cilitation of stellar flux tube model s is akin to that de- 
scribed bv lUlmschneider et al.l (|1996[ ). In this approach, in- 
formation is needed about velocities of turbulent motions 
in the overshooting layer near the stellar surface, where the 
squeezing of the magnetic flux tube is assumed to occur. 
Steffen's numerical calculations show that the rms velocities 
decrease toward the solar surface and reach a plateau in the 
overshooting layer. Between T5000 = 1 and 10 -4 , he finds 
values of w rms = 1.4 km s -1 , which are essentially indepen- 
dent of height. Ulmschncidcr & Musiclak (1998) adopted 
for the Sun a variety of observed rms velocity amplitudes 
u t in the range 0.9 < u t < 1.9 km s _1 , and showed the de- 
pendence of the computed fluxes on this velocity. For stars, 
these velocitie s cannot be determ i ned f rom observations; 
thus, we follow lUlmschneider et al. ( 2001) by assuming that 
the rms velocity fluctuations at the squeezing points are 
given by u t — «cMax/2. The values of u t and DcMax (see 
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Fig. 2) are evaluated from stellar convection zone models 
based on the adopted mixing length parameter. The nu- 
merical factor 2 used in our calculations ensures that in 
our approach the considered convective velocities are al- 
ways lower than the local speed of sound. 

After specifying the rms velocities that are responsi- 
ble for the wave generation, we also determine the height 
in stellar atmospheres, where the most efficient squeezing 
of magnetic flux tubes t akes place. We are guided by ear - 
lier studies performed bv lUlmschneider fc Musielakl |l998), 
who pointed out that shifting the height of the excitation 
point did not change much the resulting wave energy fluxes 
for the Sun. Based on these results, we take the squeezing 
point to be located at optical depth T5000 = 1 for all consid- 
ered stars. This depth is commonly taken as the zero height 
level in stellar atmosphere computations. 

For the computation of stellar magnetic flux tu bes we 
consi der the "thin flux tube approximation" (e.g., ISpruitl 
|1981[). Two dimension al modeling of magnetic flux tubes 
bv lrlasan et al.l (|2003f ) indicates that this approach renders 
reasonable results if the models do not extend beyond a 
small number of scale heights above the stellar surface; note 
that this approximation is fully compatible with our study. 
In the following, we augment th e commonly used concept o f 
solar magnetic flux tube s (e.g.,IStenflol|1978l ISolanki|[l993l) 
to stellar flux tubes; see ISolankil (|1996( ) for further discus- 
sion. Stellar magnetic flux tubes at the stellar surface are 
assumed to have diameters of roughly equal to the local 
pressure scale height like in case of the Sun. 

For the Sun, the magnetic flux tubes within the solar 
phot osphere have fi eld strengths on the order B = 1500 G 
(e.g.. lSolankl|1993D. Taking y e = 1 .17 • 10 5 dyn cm" 2 from 
model C of IVernazza et all (|l98lh at the height z = 
where T5000 = 1 yields an equipartition field strength 



of B, 



1716 G. This corresponds to a 



ratio B/B cq = 0.875, which may or may not be typical for 
stellar flux t ubes. Since this ratio is li kely to vary even on 
the Sun (e.g.. lSchriiver fc ZwaarfeOOOD . it is appropriate to 
consider a range of that ratio for the sake of comprehen- 
siveness. Therefore, we take B/B eq = 0.75, 0.85, and 0.95 
for our stellar flux tube models. This allows us to deduce 
an appropriate set of energy fluxes for upward propagat- 
ing longitudinal tube waves for specified values of a = 1.5, 
1.8, and 2.0, resulting in a total of 9 models per theoretical 
target star. 



3. Results and discussion 

3.1. Computation of wave energy fluxes 

The instantaneous and time-averaged tube wave energy 
fluxes are computed by employing a modified time- 
dependent wave cod e based on an earlier version by 
Hcrb old et al.l (|l985l) . The external turbulent motion is 
translated into internal pressure and velocity fluctuations. 
The time-averaged wave energy fluxes at the squeezing 
point are computed for each combination of effective tem- 
perature T c g, surface gravity log <?, 77 = B/B cq and mixing- 
length parameter a for flux tubes embedded in the atmo- 
spheres of the main-sequence stars as considered. 

For the evaluation of reliable time-averaged wave energy 
fluxes and due to the spiky nature of the instantaneous wave 
energy fluxes, the wave propagation code has to be run 
over a time of about 35Pd with Pjj as the wave period of 
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Fig. 4. Different types of energy fluxes (see text) for the 
set of theoretical main-sequence stars for a — 2. Results 
are given for B/B cq = 0.75 (top), 0.85 (middle), and 0.95 
(bottom). 



the Defouw cut-off frequency (|Defouwlll976h . The generated 
waves include both propagating and non-propagating wave 
energy fluxes. We apply a high pass filter at the Defouw cut- 
off frequency, wb, on the velocity and pressure fluctuations 
inside the flux tubes to separate the propagating waves and 
to compute the time-averaged upward propagating wave 
components. The position of the filter is indicated by the 
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Table 3. Wave energy flux" for different values of 77 based on a — 2 



Sp.Type T cB F(tot) -F(prop) F(prop+up) F(tot) F(prop) F(prop+up) F(tot) F(prop) F(prop+up) 
(K) ... ... 



r) = 0.75 r) = 0.85 rj = 0.95 



F2 V 


6909 










5 


18E8 


4 


37E8 


6.26E8 


1 


17E8 


1.01E8 


1.41E8 


F5 V 


6528 


8.23E8 


6 


97E8 


1.12E9 


5 


04E8 


4 


13E8 


6.07E8 


1 


27E8 


1.08E8 


1.47E8 


F8 V 


6160 


6.76E8 


6 


02E8 


8.51E8 


4 


10E8 


3 


59E8 


4.92E8 


1 


09E8 


9.48E7 


1.34E8 


GO V 


5943 


5.45E8 


4 


71E8 


6.77E8 


3 


22E8 


2 


76E8 


3.96E8 


9 


75E7 


8.11E7 


1.11E8 


G2 V 


5811 


4.91E8 


4 


26E8 


6.35E8 


2 


96E8 


2 


53E8 


3.65E8 


8 


70E7 


7.56E7 


1.06E8 


G5 V 


5657 


4.75E8 


3 


93E8 


5.81E8 


2 


73E8 


2 


32E8 


3.35E8 


8 


45E7 


7.22E7 


1.01E8 


G8 V 


5486 


3.86E8 


3 


29E8 


4.79E8 


2 


38E8 


2 


02E8 


2.86E8 


6 


76E7 


5.87E7 


8.05E7 


K0 V 


5282 


2.95E8 


2 


30E8 


3.39E8 


1 


66E8 


1 


39E8 


1.95E8 


5 


47E7 


4.68E7 


6.39E7 


K2 V 


5055 


2.52E8 


2 


08E8 


2.98E8 


1 


55E8 


1 


26E8 


1.86E8 


4 


86E7 


3.90E7 


5.81E7 


K5 V 


4487 


1.29E8 


8 


93E7 


1.44E8 


8 


08E7 


5 


93E7 


9.19E7 


2 


53E7 


1.83E7 


2.86E7 


K8 V 


4006 


5.55E6 


3 


36E6 


6.79E6 


3 


74E6 


2 


14E6 


4.37E6 


1 


25E6 


6.95E5 


1.57E6 


MO V 


3850 


2.29E6 


1 


36E6 


3.32E6 


1 


57E6 


8 


23E5 


2.11E6 


5 


37E5 


2.86E5 


7.62E5 



"The unit of F(tot), -F(prop), and F(prop+up) is erg cm 2 s 1 . Note that F(prop+up) = -Fltw- 




5 
3.85 



3.8 3.75 3.7 3.65 
log T eff (K) 



3.6 



3.55 



Fig. 5. Upward propagating energy fluxes, also referred to 
as FltW: for the set of theoretical main-sequence stars for 
a = 2 and different values of 77 = B/B eq . 



low frequency cut-off depicted in Fig. 3, which features a 
comparison of the power spectra for different types of stars. 
Note that the power spectra reach their maxima at roughly 
2ujd and decrease toward higher frequencies. 

Figure 4 and 5 as well as Table 2 and 3 show the com- 
puted wave energy of the magnetic flux tube. Actually, 
Table 3 for the case of a — 2 also distinguishes between 
the total energy flux, propa gating energy flux and upw ard 
propagating energy flux; see lUlmschneider et al.l (|200ll) for 
detailed definitions. The wave energy flux Fltw identified 
as upward propagating energy flux shows a characteristic 
behaviour as function of the model parameters, notably the 
stellar spectral type (or T e s). It is found that the value of 
Fltw decreases by about a factor of 1.7 between GO V and 
K0 V stars, and drops by almost two orders of magnitude 
between K0 V and M0 V stars. This is a direct consequence 
of the action of stellar convection, noting that in relatively 
hot main-sequence stars there is an enhanced efficiency in 



the creation of turbulence, reflected in higher values of the 
rms velocity Ut (see Fig. 2). 

Additionally, we also found that the wave energy flux 
Fltw shows significant dependence on the magnetic field 
strengths inside of the tubes, given as 77 = B/B eq as well 
as some minor dependence on the mixing-length param- 
eter a. In fact, the wave energy fluxes are found to in- 
crease with decreasing magnetic field strengths inside the 
flux tube, a phenomenon that can be explained by the de- 
creasing stiffness of the magnetic tube; see Fig. 5 for de- 
tails. This behaviour was previously also pointed out by 
lUlmschneider et alj (|2001l ). 

Moreover, it is found that the wave energy flux for up- 
ward propagating waves is also somewhat sensitive to the 
choice of a concerning the stellar convective zone. Note that 
an increase of the mixing length from a = 1.5 to 2.0 am- 
plifies the simulated longitudinal wave energy fluxes by a 
factor of approximately 1.45, which is most notable for rel- 
atively hot main-sequence stars. This can be explained by 
the fact that for higher values of a, the convective zones 
are more efficient concerning the generation of turbulence, 
thus resulting in thus higher values of the rms velocity Ut 
(see Fig. 2). 



3.2. Derivation of fitting formulas 

The amount of data for Fltw Ior stars of different spec- 
tral types and the various choices for r\ and the mixing- 
length parameter a are a strong motivation for the deriva- 
tion of a fitting formulas. Our main focus are stars be- 
tween T c ff = 6000 K (i.e., spectral type F9.5 V) and stars 
of spectral- type mid-K. Stars below mid-K are increas- 
ingly dominated by flare act ivity or other processes as, e.g ., 
magnetic reconnection (e.g.. lNarain fc Ulmschneid cr 1996), 
which makes the existence of an accurate fitting formula 
less urgent. Nevertheless, we extended our fitting formula 
toward stars of spectral type M0 V (with T e ff = 3850 K) 
while reproducing the wave energy flux of M0 V stars with 
a precision of 10~ 3 . 

However, no attempt has been made to accurately re- 
produce the knee in the wave energy flux function near F8 V 
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Table 4. Data for the parameter e 



a 


r) = 0.75 


r) = 0.85 


77 = 0.95 


1.5 


5.24 


5.30 


4.95 


1.8 


5.51 


5.11 


4.87 


2.0 


5.67 


5.31 


4.90 


Data for the parameter /j, 


Q 


77 = 0.75 


■q = 0.85 


77 = 0.95 


1.5 


4.78E-3 


4.69E-3 


4.45E-3 


1.8 


4.95E-3 


5.04E-3 


4.55E-3 


2.0 


4.47E-3 


4.53E-3 


4.49E-3 



as this would have led to an additional complication in our 
fitting formula without obvious merits. This approach is 
motivated by the finding that energy dissipation by longi- 
tudinal tube waves appears to be less important in dwarfs 
stars of spectal type mid -K to M compared to, e.g., flare 
heating as pointed out bv lFawzv et al.l (|2002j ) in their com- 
parison between empirical and theoretical radiative chro- 
mospheric emission losses. 

Our formula is given as follows: 



LTW 



= 2.55 • 10 s Tt 



f Z(a, rj) erg cm 2 s 



(1) 



with T re i = T cff /T effi0 noting that T cffj © 
addition, Z(a,r)) is given as 



Z(a,rj) 



, 1.25 



V0.75/ 



5777 K. In 



(2) 



with r\ = B/B cq and 7 = —4.5 — 30- \rj — 0.85 1 . Furthermore, 
Tdif is defined as Min(0, T c g — 4500 K), which means that 
Tdif is zero for T e g > 4500 K and takes a negative value 
otherwise. Note that the parameters e and /i of Eq. (1) 
weakly depend on a and 77; see Table 4 and 5 for detailed 
information. If a reduced level of accuracy is permitted, it 
might be appropriate to use average values for a and /z, 
given as a — 5.21 and (1 — 4.66 • 10~ 3 , respectively. 

Our formula for the upward propagating longitudinal 
wave energy flux has been subjected to thorough testing 
for stars of spectral type F9.5 V, G2 V, G5 V, G8 V, K0 V, 
K2 V, and K5 V. The simulated data for the F9.5 V star 
were obtained via logarithmic interpolation between the 
data for the F8 V and GO V stars. Detailed information 
on the tests is given in Appendix A. 



4. Summary and conclusions 

We studied the generation of longitudinal waves in stellar 
magnetic flux tubes of theoretical main-sequence stars. Our 
results are commensurate with those obtained from pre- 
vious studies, especially the work by lUlmschneider et al.l 
(2001). Our investigations show that through nonlinear 
time-dependent responses of stellar magnetic flux tubes to 
continuous and impulsive external turbulent pressure fluc- 
tuations, longitudinal tube waves were effectively produced 
via dipole emission. Furthermore, the shapes of the com- 
puted power spectra were found to be similar for stars of 
different effective temperature. Moreover, the longitudinal 



wave energy fluxes are found to increase with higher effec- 
tive temperature, i.e., stars of earlier spectral types. 

As part of our study, we investigated the role of the mag- 
netic field strength inside of the tube B as well as that of 
the adopted convective model characterized by the mixing- 
length parameter a concerning the generated wave energy 
flux. We found that the computed wave energy flux strongly 
depends on the strength of the m agnetic field as already dis- 
cussed in previous studies (e.g ., lUlmschneider fc Musielald 
119981: lUlmschneider et aLll2001h . For a given spectral type, 
the flux is considerably higher in tubes with a field strength 
of B/B cq = 0.75 compared to B/B cq = 0.95, although the 
difference as a function of spectral type is not as l arge as 
previously pointed out bv lUlmschneider et al.l (|2001l) owing 
to the differences in stellar surface gravity for the different 
types of stars. Note that the difference concerning the mag- 
netic field strength of B/B cq = 0.75 and 0.95 is found to be 
a factor of 7.6, 6.1, 5.3, and 4.4 (for a = 2) for stars of spec- 
tral type F5 V, GO V, K0 V, and M0 V, respectively. This 
difference exhibits a noticeable, albeit little dependence on 
the mixing-length parameter a for stars hotter than G2 V, 
rooted in the behaviour of the adopted root mean square 
velocity at the squeezing point of the tube. 

Another aspect of our study was to consider a limited 
range of the mixing-length parameter a, which are 1.5, 1.8, 
and 2.0. It was found that an increase of the mixing length 
from a = 1.5 to 2.0 enhances the computed energy fluxes by 
a factor of about 1.45, corresponding to a proportionality 
of a 1 ' 25 . This relationship can be compared with previous 
findings for acoustic e nergy gener ation, which show a de- 
p ende nce such as cn 8 (lBohnlll984ft or a 3 ' 8 in updated mod- 



els bv lMusielak et al.l (|1994 ). The relatively weak influence 
of a on the amount of upward propagating wave energy flux 
-Fltw is apparently due to the fact that the latter is largely 
controlled by magnetic processes, as also evidenced by the 
strong influence of 77 up to 77~ 7 ' 5 , rather than controlled by 
convective processes, even though the latter are key for the 
excitement of the tubes subsequently resulting in magnetic 
wave generation. 

Note that the strong dependence of the generated wave 
energy fluxes on the stellar and magnetic parameters is in 
general agreement with the findings from previous stud- 
ies, although some noticeable differences are attained. The 
wave energy flux -Fltw decreases by about a factor of 1.7 
between the GO V and K0 V stars, and drops by almost two 
orders of magnitude between the K0 V and M0 V stars. For 
a = 2 at a fixed value of T g = 5000 K, Ulms chneider et all 
(2001) deduced a fitting formula for the behaviour of Fltw 
as function of B/B eq , which is found to be commensurate 
with the results obtained in our current study. On the other 
hand, the fitting formula for the wave energy flux given in 
our paper is more general than any of the previously de- 
duced formulas as it is applicable to a large range of stellar 
effective temperatures and furthermore allows insight into 
the role of the governing magnetic and convective parame- 
ters concerning the amount of generated upward propagat- 
ing wave energy flux. Therefore, it is of interest to future 
solar and stellar physics studies as it allows flexibility both 
concerning the mixing-length parameter a and the mag- 
netic parameter 77 = B/B cq . Note that even for the Sun, 
77 is expected to exhibit considerable spatial and temporal 
fluctuations across the surface as implie d by previous stud- 
i es of solar physics research; see, e.g., ISchriiver fc Zwaaiil 
(2000) for background information. 
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Table A.l. Accuracy of the fitting formula 



References 



a 


Metric" 


r) = 0.75 


7} = 0.85 


T) = 0.95 


1.5 


devi 


5.2 


6.1 


5.6 


1.5 


dev2 


5.8 


7.5 


6.2 


1.5 




9.5 


15.0 


9.1 


1.8 


devi 


6.1 


3.4 


2.9 


1.8 


dev2 


7.1 


4.0 


3.6 


1.8 




11.7 


6.5 


6.5 


2.0 


devi 


6.9 


7.3 


5.8 


2.0 


dev2 


7.3 


7.9 


7.3 


2.0 




12.4 


10.1 


15.5 



"Data for devi, dev2 and A max are given in percent. 



Appendix A: tests for the fitting formula 

For testing our fitting formula (see Eq. 1) we considered two 
different metrices, i.e., the linear and the quadratic metric, 
which allow us to assess the accuracy of the formula. The 
linear metric is given as 



1 N 

deVl = JjYl I F ltw - -Pltw I , 



(A.l) 



whereas the quadratic metric (also referred to as rms met- 
ric) is given as 



dev 2 = 



\ 



1 N 

i = l 



i 7 ! 



LTWl 



(A.2) 



Here -Fltw refers to the wave energy flux of the detailed 
models, whereas -F£ TW refers to the wave energy flux ob- 
tained by the fitting formula (or vice versa). N denotes the 
number of stars per test series. 

The test results have been deduced in a separate manner 
for the different values of 77 and a; see Table A.l for details. 
It is found that the average deviation, regardless of the 
selected metric, are typically considerably better than 10%, 
and for some of the test series, the average deviation is 
found to be better than 5%. We also checked the maximal 
deviation for individual stars, denoted as A max , for a given 
test series. Our results indicate that the maximal deviation 
never exceeds 15%. Finally, we also calculated the mean 
deviation between the model data and the data given by the 
formula for the entire set of theoretical main-sequence stars, 
comprised of 63 models. We found that for the entire set 
of model stars, the linear metric yields a mean deviation of 
5.5%, whereas the quadratic metric yields a mean deviation 
of 6.5%, a strong testimony of the quality of our fitting 
formula for the overall range of solar-type stars. 



Acknowledgements. This work has been supported by the Faculty of 
Engineering and Computer Sciences, Izmir University of Economics 
(D. E. F.) and the Department of Physics, University of Texas at 
Arlington (M. C). The authors also appreciate previous comments 
by P. Ulmschneider and Z. E. Musielak. 



Bohn, H. U. 1984, A&A, 136, 338 

Buchholz, B., Ulmschneider, P., & Cuntz, M. 1998, ApJ, 494, 700 
Castelli, F., & Kurucz, R. L. 2004, in Modelling of Stellar Atmospheres 

ed. N. E. Piskunov, W. W. Weiss, & D. F. Gray (San Francisco: 

ASP), IAU Symp. 210, CD-ROM, Poster 20 
Cattaneo, F., Brummell, N. H., Toomre, J., Malagoli, A., & Hulburt, 

N. E. 1991, ApJ, 370, 282 
Cuntz, M., Ulmschneider, P., & Musielak Z. E. 1998, ApJ, 493, L117 
Cuntz, M., Rammacher, W., Ulmschneider, P., Musielak, Z. E., & 

Saar, S. H. 1999, ApJ, 522, 1053 
Cuntz, M., Rammacher, W., & Musielak, Z. E. 2007, ApJ, 657, L57 
Cuntz, M., Guinan, E. F., & Kurucz, R. L. 2010, in Solar and Stellar 

Variability: Impact on Earth and Planets, ed. A. G. Kosovichev, A. 

H. Andrei, & J. -P. Rozelot (Cambridge: Cambridge Univ. Press), 

IAU Symp., 264, 419 
Defouw, R. J. 1976, ApJ, 209, 266 
Fawzy, D. E. 2010, MNRAS, in press 

Fawzy, D. E., Ulmschneider, P., & Cuntz, M. 1998, A&A, 336, 1029 
Fawzy, D., Ulmschneider, P., Stcjpieh, K., Musielak, Z. E., & 

Rammacher, W. 2002, A&A, 386, 983 
Gray, D. F. 2005, The Observation and Analysis of Stellar 

Photospheres, 2nd edn. (Cambridge: Cambridge Univ. Press) 
Giidel, M. 2007, Liv. Rev. Sol. Phys., 4, 3 
Guinan, E. F., Ribas, I., & Harper, G. M. 2003, ApJ, 594, 561 
Hasan, S. S., Kalkofen, W., van Ballegooijen, A. A., & Ulmschneider, 

P. 2003, ApJ, 585, 1138 
Herbold, G., Ulmschneider, P., Spruit, H. C, & Rosner, R. 1985, 

A&A, 145, 157 
Komm, R., Mattig, W., & Nesis, A. 1991, A&A, 243, 251 
Kurucz, R. L. 2005, Mem. S. A. It., 8, 14 

Lammer, H., Selsis, F., Ribas, I., Guinan, E. F., Bauer, S. J., & Weiss, 

W. W. 2003, ApJ, 598, L121 
Muller, R. 1985, Sol. Phys., 100, 237 

Muller, R., Roudier, Th., Vigneau, J., & Auffret, H. 1994, A&A, 283, 
232 

Musielak, Z. E., Rosner, R., & Ulmschneider, P. 1989, ApJ, 337, 470 
Musielak, Z. E., Rosner, R., Gail, H. P., & Ulmschneider, P. 1995, 
ApJ, 448, 865 

Musielak, Z. E., Rosner, R., Stein, R. F., & Ulmschneider, P. 1994, 
ApJ, 423, 474 

Musielak, Z. E., Rosner, R., & Ulmschneider, P. 2000, ApJ, 541, 410 
Narain, U., & Ulmschneider, P. 1990, Space Sci. Rev., 54, 377 
Narain, U., & Ulmschneider, P. 1996, Space Sci. Rev., 75, 453 
Nesis, A., Hanslmeier, A., Hammer, R. , Komm, R. , Mattig, W., & 

Staiger, J. 1993, A&A, 279, 599 
Nordlund, A., & Dravins, D. 1990, A&A, 228, 155 
Nordlund, A., & Stein, R. F. 1991, in Stellar Atmospheres: Beyond 

Classical Models, ed. L. Crivellari & I. Hubeny (Dordrecht: 

Kluwer), 263 

Nordlund, A., Spruit, H. C, Ludwig, H.-C, & Trampedach, R. 1997, 
A&A, 328, 229 

Saar, S. H. 1994, in Cool Stars, Stellar Systems, and the Sun 8, ed. 

J.-P. Caillault (San Francisco: ASP Conf. Ser.), 64, 319 
Schrijver, C. J. 1996, in Stellar Surface Structure, ed. K. G. 

Strassmeier & J. L. Linsky (Dordrecht: Kluwer), IAU Symp. 176, 

1 

Schrijver, C. J., & Zwaan, C. 2000, Solar and Stellar Magnetic Activity 

(Cambridge: Cambridge Univ. Press) 
Schroder, K.-P., & Eggleton, P. P. 1996, Rev. Mod. Astr., 9, 221 
Solanki, S. K. 1993, Space Sci. Rev., 63, 1 

Solanki, S. K. 1996, Stellar Surface Structure, ed. K. G. Strassmeier 

& J. L. Linsky (Dordrecht: Kluwer), IAU Symp. 176, 201 
Spruit, H. C. 1981, A&A, 98, 155 

Steffen, M. 1993, Habilitation Thesis, Univ. Kiel, Germany 

Stein, R. F., Nordlund, A., Georgobiani, D., Benson, D., & 
Schafenberger, W. 2009a, in Solar-Stellar Dynamos as Revealed 
by Helio- and Asteroseismology, ed. M. Dikpati, I. Gonzalez- 
Hernandez, T. Arentoft, & F. Hill (S an Francisco: ASP ), GONG 
2008 / SOHO XXI Meeting, in press llarXiv:661TT 0472 I 

Stein, R. F., Georgobiani, D., Schafenberger, W., Nordlund, A., & 
Benson, D. 2009b, in Cool Stars, Stellar Systems, and the Sun 15, 
ed. E. Stempels (Melville: AIP Conf. Proc), 1094, 764 

Stenflo, J. O. 1978, Rep. Progr. Phys., 75, 3 

Trampedach, R., Christensen-Dalsgaard, J., Nordlund, A., Stein, 
R,. F. 1997, in Solar Convection and Oscillations and their 
Relationship, ed. F. P. Pijpers, J. Christensen-Dalsgaard, & C. S. 



Diaa E. Fawzy and M. Cuntz: Longitudinal flux tube 



Rosenthal (Dordrecht: Kluwer), 73 
Ulmschneider, P., & Musielak, Z. E. 1998, A&A, 338, 311 
Ulmschneider, P., Theurer, J., & Musielak, Z. E. 1996, A&A, 315, 212 
Ulmschneider, P., Theurer, J., Musielak, Z. E., & Kurucz, R. 1999, 

A&A, 347, 243 

Ulmschneider, P., Musielak, Z. E., & Fawzy, D. E. 2001, A&A, 374, 
662 

Vernazza, J. E., Avrett, E. H., & Loeser, R. 1981, ApJS, 45, 635 



